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We report on experimental results for the conductance of metallic single-electron transistors as 
a function of temperature, gate voltage and dimensionless conductance. In contrast to previous 
experiments our transistor layout allows for a direct measurement of the parallel conductance and 
no ad hoc assumptions on the symmetry of the transistors are necessary. Thus we can make a 
comparison between our data and theoretical predictions without any adjustable parameter. Even 
for rather weakly conducting transistors significant deviations from the perturbative results are 
noted. On the other hand, path integral Monte Carlo calculations show remarkable agreement with 
experiments for the whole range of temperatures and conductances. 

PACS numbers: 73.23.Hk, 85.35.Gv, 02.70.Ss 



I. INTRODUCTION 

The usual single-electron transistor (SET) layoutQ! 
consists of two ultrasmall tunnel junctions with conduc- 
tances G s ,d and capacitances C s ,d, respectively. The 
junctions in series are biased by a voltage C/ S( j and the 
island between the tunneling barriers is coupled to a gate 
voltage U g via a capacitance C g . 



as the basic device to 
Its transport properties 



This setup has been us 
study the Coulomb blockad 
are governed by two dimensionless parameters. Firstly, 
the dimensionless inverse temperature (3E C relates (3 — 
(ksT)^ 1 to the charging energy E c = e 2 /(2Cs) with 
Cs = C s + Cd + C g being the total island capaci- 
tance. This parameter determines how far the electro- 
static blockade of the source-drain current is lifted by 
thermal excitations. Secondly, the dimensionless parallel 
conductance g = (G s + Gd)/Gx, Gk — e 2 /h being the 
conductance quantum, measures how much the quanti- 
zation of the island charge is smeared by quantum fluc- 
tuations. The charge on the island can be controlled 
by means of the gate voltage U g . The linear-response 
conductance G between source and drain is a periodic 
function of U g taking on values between G m i n and G max - 
Since the control of a device by Coulomb blockade ex- 
ploits a large difference between G m in and G max a de- 
tailed understanding of the washout of Coulomb blockade 
effects by thermal and quantum fluctuations is crucial for 
the optimal design of fast and reliable Coulomb blockade 
devices. 

Theoretical work has mainly focused on the limits of 
small g or f3E c , respectively. In the-JKeakly conducting 
regime, g <C 1, perturbation theoryotffl should be suffi- 
cient to describe the experimental data. For (3E C <C 1 it 
is more suitable to formulate the problem in terms rjafi 
path integral which may be evaluated semiclassicallyOEH 
Only recently the use of path integral Monte Carlo 
(PIMC) techniques was proposed to calculate the con- 



ductance of the SET over tbe whole range of experimen- 
tally accessible parametersEl Especially the regime of 
large parallel conductance g is of interest for technologi- 
cal applicationsOHLil 

This recent work has pointed to a shortcoming of previ- 
ous experiments. The first step of a comparison between 
theory and experiment is the determination of the pa- 
rameters g and E c . For the SET layout described above 
one can obtain only the series conductance Gs of the 
two junctions which does not suffice to calculate the di- 
mensionles s parallel co inductance g without further as- 



Joyez et al 



assumed that their SETs were 



sumption. 

symmetric, i. e. built up of two identical tunnel barri- 
ers which turned out to be a good approximation for the 
first three samples measured in the experiment. However, 
comparing the data for the high condu ctance sample 



Rcf. [j]to their PIMC simulation results, Soppcrt et al 



found strong indications that the tunnel barriers could 
be asymmetric. 

Since an unambiguous determination of the relevant 
parameters is a prerequisite of a thorough comparison 
between theoretical and experimental results we have de- 
veloped a transistor design with four tunnel junctions 
connected to the island (see Fig. |l|). This arrangement 
allows the determination of the individual resistances of 
each junction and therefore the direct measurement of 
the tunneling strength parameter g. It is operated as 
SET by connecting two of the junctions in parallel to 
form source and drain, respectively. 

In Sec. H we present experimental details about the 
sample fabrication, the determination of the parameters 
and the conductance measurements. In Sec. Ill we give 
the path integral formulation which is used for the imag- 
inary time quantum Monte Carlo simulation. We sum- 
marize how the conductance can be calculated from the 
simulation data by use of the singular value decomposi- 



tion (SVD) analytical continuation scheme. In Sec. IV 



we compare experimental and theoretical results. At low 
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(a) 




FIG. 1: a) Circuit diagram for the four junction layout when 
operated as SET. b) SEM picture of a four junction SET. 
The layout contains two symmetrically arranged gate fingers 
to avoid asymmetrical proximity effects in the lithography 
process. In the experiments both gates are connected in par- 
allel. 



conductance the perturbation theory is sufficient to de- 
scribe the temperature dependence of G max , but for G m ; n 
significant deviations are observed. At high conductance, 
we make a direct comparison of the experimental data 
with path integral Monte Carlo results. 



II. EXPERIMENT 

Seven samples with varying tunneling strengths were 
investigated. The samples were fabricated from alu- 
minum by standard e-beam lithography in combination 
with two-angle shadow evaporation. The evaporation is 
performed by electrical heating a tungsten wire which 
holds a drop of aluminum. Tunnel barriers of different 
strengths could be produced by a variation of the oxygen 
pressure applied to the evaporation chamber between the 
two evaporation steps. 

We used two different layouts: firstly, the usual SET 
design with two junctions, forming a small island in be- 
tween, and a straight gate finger from the side pointing 
towards the island, secondly, the design which is shown 



on the SEM picture in Fig. |T|. In this design one can 
determine the individual tunnel resistances by measur- 
ing the current in response to an applied voltage bias 
across different combinations of the four tunnel junc- 
tions. The measurements presented in this article are 
performed on four standard SETs (samples I-IV), one 
four-contact sample (sample VII) and two samples which 
also have the four-contact layout, but turned out to have 
only three working tunnel junctions (samples V and VI). 
Nevertheless, this is sufficient to determine all individ- 
ual tunnel resistances. The sample parameters arc given 
in Tabic ||. The measurements are performed in a top- 
loading dilution refrigerator in the temperature range 
from 25 mK to 18 K. The samples are mounted within a 
well-shielded metallic cavity. All electrical wiring into the 
cavity is made of highly resistive leads (32 f2/m, diameter 
0.23 mm) which are fed through stainless steel capillaries 
with an inner diameter of 0.34 mm and a length of 1 m. 
The capillaries are wound up in a compact coil which 
is held in thermal equilibrium with the sample. These 
feedthroughs form resistive co-axial cables. They pro- 
vide a damping exceeding 200 dB in the frequency range 
from 20 GHz to 6 THz, as calculated by classical elec- 
trodynamics taking the skin effect into accounto. The 
validity of the used formulae was checked experimentally 
in the frequency range up to 20 GHz using a spectrum 
analyzer. Additional rf-filtering is performed at room 
temperature at the entrance to the cryostat. 



To measure the conductance of the transistors they are 
biased with a voltage and the resulting current is mea- 
sured with an operational amplifier at the top of the cryo- 
stat. The resolution of the current measurement is bet- 
ter than 100 fA. To gain resolution and to circumvent the 
1 /f- increase of noise at low frequencies an AC component 
of w 10 Hz is added to the biasing DC voltage simultane- 
ously and the resulting AC component of the current is 
measured with a lock-in amplifier. The DC measurement 
is used to ensure that the measurement is performed at 
vanishing bias and stays in the linear-response regime. 



To determine the conductance Gs of two tunnel junc- 
tions in series, we measure the I s dU s d characteristic up to 
a bias voltage of maximal ±20 mV. We define the asymp- 
totic slope at large bias voltages as G^. The four-junction 
design allows six different configurations of two tunnel 
junctions connected in series. From the corresponding set 
Gs^, i = 1, 6, the individual tunnel junction conduc- 
tances can be derived with simple algebra. Their values 
are given in table |. Thus this layout enables us to actu- 
ally measure the coupling strength parameter g directly. 
The SET investigated in the following experiments was 
formed by using all four tunnel junctions, where source 
and drain were made by connecting two tunnel junctions 
in parallel, respectively (see Fig. ^). 
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VI 


59.4 
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18.0 


22.7 




104.2 


56.1 


2.7 


1.85 


VII 


23.0 


497 


19.0 


20.3 


16.4 


31.7 


23.8 


1.87 


4.75 



TABLE I: Parameters of the samples I- VII. G c i denotes the high temperature conductance of the SET. The gate capacitance 
C g is extracted from the period of the Coulomb oscillations within an accuracy of about 1%. Samples V-VII consist of at 
least three tunnel junctions. The individual conductances G a , a > and Gd,d' for these samples were evaluated by simple algebra 
from values of Gs, which had been measured at different pairs of tunnel junctions as described in the text. For the 2-junction 
samples I-IV, the given value of g equals (G s + Gd) /(2Gk), an expression valid for symmetric SETs only. For the asymmetric 
sample I, the value derived from a comparison with the second-order perturbation theory is given in brackets (see Fig. 



III. THEORY 



with the partition function 



A. Path integral formulation 

The linear DC conductance G represents a transport 
coefficient which can be expressed by correlation func- 
tions of the system using a Kubo formulaM Defining 
the current through the SET as the average of the current 
through the source and drain junctions I = (J s + Id) /2 
the conductance G may be connected to the spectrum of 
the current autocorrelation function F(t) — (Z(t)I(O)), 
i.e. 



G=^F(uj = 0) 



(1) 



with F(cu) denoting the Fourier transform of F(t). The 
calculation of the cupsent correlator may be done in the 
phase representationta, i.e. in terms of the phase variable 
ip which is conjugate to the charge q on the island. For 
imaginary time r one getsB 

F(t) = 47rG cl a(r)C(r) (2) 

with the classical high temperature conductance G c \ and 



Z7T 

a(u>) 



dtu a{oj) e 



— oo 



27T 1 - e~ h ^ ' 



(3a) 
(3b) 



The cosine correlation function C(t) = 
(cos [(p(r) — <p(0)]) has the formally exact path in- 
tegral representatiorfl 

^ <p(0ti)=2wk 

fe=_0 ° v (o)=o 

x cos [<p(r) - tp(0)] (4) 



Z 



i., j 



rS[ip] e 2irikn g 



(5) 



and the dimensionless gate voltage n g — (U g C g )/e. The 
Euclidian action S[(p] — Sc[f] + S'tIv] splits into the 
Coulomb action 



Sr 



o ±E C 



(6) 



which describes the charging of the island and the tun- 
neling action 



0H ph 

S T [<p] = 2.g J dr J dr 1 a(r - r') sin 2 
o o 



tp(r) - ip(t') 



(7) 



that expresses the influence of tunneling processes on the 
dynamics of the phase variable <p. 

Equation (Q) W&y serve as the starting point for analyt- 
ical calculationa§EI or numerical worklld. In the latter case 
the interval [0, f3h] is divided into N Trotter slices and 
the multidimensiQiial integral is calculated using Monte 
Carlo methodsESLj. Since the action S[<p] is real, the 
Metropolis algorithm can be applied for an importance 
sampling of the configurations {tpi = (p(ri)\i — 0, ...,iV} 
and the winding number k. 

We did Monte Carlo simulations for fixed tunneling 
strength g = 4.75 over a range of inverse temperatures 
j3E c e [0.5,21.0] which will be compared to our exper- 
imental findings in Section PV. For each temperature 
the system was equilibrated during several million Monte 
Carlo steps. Measurements were then carried out for an- 
other 5-32 million sweeps depending on the value of (3E C . 
Especially for low temperatures it is necessary to increase 
the number of measurements to get reliable statistics of 
the data. We ensured that the error of the correlation 
function is always less than 3%. Over the whole range 
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of temperatures we chose N = 200 Trotter slices fulfill- 
ing t he , c onvergence criterion N > 5f3E c used in earlier 
workllllla. The conductance was calculated for 200 val- 
ues of the dimensionless gate voltage spanning the range 
n g £ [0,0.5]. 

Additionally we examined single electron transistors 
with different parallel conductances g £ [2.0, 15.0] for 
fixed temperature 1/((3E C ) — 0.05. Here we chose N = 
250 Trotter slices and did up to 23 million measurement 
sweeps. We found that the convergence was slower for 
small values of the dimensionless conductance. 



B. Analytic continuation 

Having calculated the cosine correlation function C (r) 
using the PIMC method we still have to solve Eqs. (0) 
and ^ for the linear conductance ij. With the posi- 
tive and symmetric spectral function!!!] A(u>) = C(u))(l — 



with the coefficients 



)/lo one can write the Fourier transformation of 



C(r) in imaginary time as 

C = KA 

where the integral operator K is given by 



iKA)(r) = - 



uj cosh 



diu 



\ 2 



r)u> 



smh(^w) 



(8) 



A(u>) (9) 



and Eqs. (0) and (|J) can be combined to give 



2n 



n cosh((3huj) 



(10) 



Before we can calculate G using Eq. ( jT0| ) it is necessary 
to invert the integral transform of Eq. (g|) to obtain A(ui). 
This problem is similar to the analytical continuation of a 
function from the imaginary to the real axisI3. It belongs 
to the class of ill-posed inverse problems. A straightfor- 
ward inversion of Eq. (||) fails, because the integral op- 
erator K is almost singular and the statistical error in 
the data C(r) is strongly amplified, making the result 
A(u>) meaningless. To deal r*fith this problem, we use 
a recently developed methodcj based, upon the singular 
value decomposition of the integral operator. 

From the SVD we get the singular system of K which 
fulfills 

{Kuj)(u) = o-jVjiu) (11a) 
{K^v 3 ){t) = a jUj (r) (lib) 

O-Q > (71 > CT 2 > ... > . (lie) 

The functions Uj(r) and Vj(uj) are called the right and 
left singular vectors, respectively. The real and positive 
a j are the singular values of K. The formal solution to 
the inverse problem (0) is given asE3 



— Pi 

j=0 J 



(12) 



tin 



dr C(r) Uj (r). 



(13) 



Since the numerical calculation of the correlation func- 
tion C(t) can only be done for a discrete set of points 
and additionally introduces a statistical error, the expan- 
sion coefficients (jl3| ) have a limited accuracy. Taking into 
account that for an ill-posed problem such as (||) the sin- 
gular values Uj vanish overexponentially with increasing 
j it becomes obvious that only the first few terms in the 
expansion ( |l2[ ) contain meaningful information whereas 
the higher terms will just corrupt the result by amplify- 
ing the noise of the data. This idea is used in the trun- 
cated SVD approachE3E2l which truncates the summation 
in ( [l2] ) by neglecting all terms for which o~j / <jq is smaller 
than the statistical error of the correlation function C (r) . 

Apart from Eq. (||) we still have supplementary infor- 
mation, namely that the spectral function A(u>) is pos- 
itive and symmetric. To £asure positiveness one may 
use a "triangular window"Lil, i.e. in the truncated sin- 
gular value decomposition one multiplies the singular 
values with a weight factor wh i ch falls of linearly from 
one to zero. Previous studiestSEd have shown that this 
implementation of positiveaess reduces the resolution of 
the method. Recent worked has shown how to use sup- 
plementary information of positiveness to enhance the 
resolution of the singular value decomposition method. 
The idea is to determine additional expansion coefficients 
which cannot be inferred from the inverse problem. They 
can be fixed by the constraints that the result shall be 
positive and the difference to the truncated SVD solution 
shall be minimal. Further details about the implementa- 
tion of the method can be found in Ref. ^0[ 

The approach used in this work is limited only by the 
statistical error in the correlation function (^J). The er- 
ror of our PIMC calculation and the subsequent ana- 
lytical continuation was estimated as follows. First we 
determined the statistical error of the Monte Carlo data. 
Then we produced an ensemble of data sets with different 
realizations of the error by adding Gaussian distributed 
random noise of the given size. The error bars shown in 
Figs. ^, ^ and ^ were then calculated from the maximum 
and minimum result of the analytical continuation. 



IV. RESULTS AND DISCUSSION 

In this section we present experimental data and the- 
oretical results for the minimum and maximum linear 
response conductance, G min and G max , respectively, as a 
function of the dimensionless parallel conductance g and 
the temperature. Before we turn to the comparison be- 
tween experimental and theoretical data we would like to 
summarize our findings about the determination of the 
experimental parameters. 

The four-junction layout allows us to determine the di- 
mensionless parallel conductance unambiguously offering 
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FIG. 2: Gmax and G m in normalized to the high temperature 
conductance G c i for sample I (•) and II (o) as a function of the 
normalized temperature. Gmax and Gmin are the maximum 
and minimum linear response conductance observed as a func- 
tion of the gate voltage U g , respectively. Lines correspond to 
the predictions of a second-order perturbative expansion in g 
(solid line: g = 1.4, dashed line: g = 1.1, dashed-dotted line: 
g = 0.8). The data are best described by g = 1.1 (•) and 
<? = 1.4(c). 



the possibility of a direct comparison between experiment 
and theory. On the other hand measurements carried out 
on the two-junction samples I and II (see Fig. |^) clearly 
demonstrate a problem also encountered in Rcf. [L0]. The 
assumption that the conductance is distributed equally 
among both tunnel junctions would lead to g\ = 0.80 
and gu = 1.39 for the samples I and II in contradiction 
to their almost identical temperature dependence of both 
G m i n and G max . Since the maximum conductance G max 
in the weak tunneling regime should ikll off proportional 
to ghi(l3E c /TT) for low temperaturescl Fig. || shows that 
the correct parallel conductances of those samples are al- 
most equal. Our four-junction samples allow us to check 
this assumption directly (see below). We found that the 
tunnel junctions are not identical although they are pro- 
duced simultaneously during shadow evaporation. This 
is not surprising as we tried to produce contacts with 
oxide barriers as thin as possible and very small junc- 
tion areas. At the borderline between functioning and 
broken contacts the unavoidable variations in the fabri- 
cation procedure become visible as large fluctuations in 
the conductances of different junctions. 

In previous experimentsBSIil several methods for the 
determination of the charging energy E c were proposed. 
As already mentioned in Ref. || the determination of the 
charging energy from the offset of the I s dU s d curve is 
not very accurate. For sample V and VII we have ana- 
lyzed the subgap resonances observed in the SET in the 
superconducting state to obtain the rcnormalizcd charg- 
ing energy as described in Rcf. [l]. |Joycz ct~aT also 
give a perturbative [tesult which connects E^ to the bare 
charging energy EjEB. Unfortunately it is only valid up 



G c , - G 




FIG. 3: Inverse conductance reduction G/(G c i — G) as a 
function of temperature. (•): Experimental data for sam- 
ple VII, measured at vanishing source-drain voltage. ( — ): 
High-temperature expansion Eq. (jl4j) with E c chosen so as to 
fit the data. 



to 0{g 2 ). In our experiments we found that this is not 
sufficient for a description of high conductance single- 
electron transistors. In our opinion the best method for 
the determination of the charging energy is a compari- 
son of the high-temperature experimental data to semi- 
classical calculations!! Apart from the full semiclassical 
expressions we tested the high-temperature expansion 



G 



ik B T 27.gC(3) 



G c i — G 



2tt 4 



(14) 



using only data with ksT ^S> gE c /2'K A (see Fig. ||). At 
a dimensionless conductance of g = 4.75 the difference 
in the charging energy between the full semiclassical re- 
sult andjEq. ( |i"4| ) was less than 1%. In contrast to ear- 
lier worku sufficient experimental data at high tempera- 
tures are available for a determination of E c . Moreover, 
Eq. ( [14] ) provides a consistency check for the dimension- 
less parallel conductance g which we have determined 
independently. 

Having determined the experimental parameters we 
can turn to the comparison with the theoretical results. 
First of all we want to compare the data measured on 
sample I, III and V-VII with perturbation theory in sec- 
ond order in the parallel conductance g. Fig. shows 
the maximum and minimum conductance of the single- 
electron transistor as a function of temperature. For sam- 
ple I and V with their moderate g = 1.1 and g = 1.40, 
respectively, good agreement can be stated for G max . At 
higher g (sample III, VI and VII) deviations of increas- 
ing size are visible. Such deviations are not surprising as 
the perturbation expansion is not justified in this param- 
eter range. For the minimum conductance G m i n we get 
deviations from the second-order result for all values of 
g investigated. At the lowest temperatures the determi- 
nation of the minimum conductance G m i„ is limited by 
the resolution of the lock-in signal which is of the order 
10 _3 Go. However the deviations are also prominent at 
1/((3E C ) = 0.05 where the lock-in signal has a sufficient 
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FIG. 4: Maximum (a) and minimum (b) conductance normalized to the high temperature conductance G c \ for samples I, V, 
VI, VII, and III with (from top to bottom) g — 1.10, 1.40, 1.85, 4.75, and 5.38 as a function of the normalized temperature. 
Together with the experimental data (•) the predictions of the perturbation theory in second order ( — ) are shown for all 
samples. For sample VII the results of the Monte Carlo calculations (□) are given additionally. For the sake of clarity the 
curves are plotted in (a) with a vertical offset increasing by 0.2 from curve to curve, in (b) the different datasets are multiplied 
by 10 2 , 10 4 , etc. 



accuracy. The discrepancy between perturbation theory 
and our data increases with g, indicating that higher- 
order corrections hawa to be included. With third-order 
perturbation theoryo deviations for G m i n are reduced 
but still visible. 

Besides renormalization group methods^ the quantum 
Monte Carlo approach is the only method which can 
cover the whole range of parameters that is accessible 
in the experiment. For sample VII (<? = 4.75), which is 
beyond the perturbative regime, we perform a detailed 
comparison. 

In Fig. H we show the gate voltage dependent conduc- 
tance for inverse temperatures ranging from 1/((3E C ) = 
0.048 to l/(j3E c ) = 2.0. We find that the experimen- 
tal Coulomb oscillations are very well described by the 
Monte Carlo calculations. Minor discrepancies at some 
temperatures can be attributed to the fact that the tem- 
peratures used in the simulations do not exactly match 
those of the experimental data. For the lowest temper- 
atures of the experiment it was not possible to get con- 
verged Monte Carlo results for the whole range of gate 
voltages in reasonable time. 

In Fig. ^ the data are analyzed in terms of the mini- 
mum and maximum conductance. Once again the accor- 
dance is remarkable. For high (3E C the maximum conduc- 
tance could not be determined by our Monte Carlo calcu- 
lations. Here a limitation of the Monte Carlo procedure 
becomes obvious. For low temperatures more terms of 
the winding number summations in Eqs. (Q) and (||) are 



relevant leading to phase cancellations due to the factors 
cxp(2irikn g ) which are especially strong at n g = 0.5, i. e. 
for the maximum conductance. Thus the convergence of 
the Monte Carlo procedure gets slower with decreasing 
temperature and the data cannot be determined as ac- 
curately. Since the analytic continuation is sensible to 
the statistical error of the data, reliable results for G max 
could not be obtained with reasonable effort and the data 
points for G max at the lowest temperatures have been 
omitted. This is not the case for the minimum conduc- 
tance G m i n as can be seen in the inset of Fig. |^. Here 
no phase cancellations occur and the experimental and 
theoretical data match nicely even on a logarithmic plot. 
We can also observe that in contrast to the perturba- 
tion theory in second order (cf. Fig. ||), the Monte Carlo 
approach gives an accurate description of G m i n at low 
temperatures. 

Finally we have examined the maximum and minimum 
conductance of single-electron transistors for varying tun- 
neling strength g at a fixed temperature l/(f3E c ) = 0.05. 
The results are show n in Fig, j?] . Besides our experimen- 
tal data the results of Joyez et al. U are shown. According 



to Ref. |10| the data of Joyez et al 
displayed at g = 10 



the earlier results of Goppert et al 



at g — 7.5 has been 
The Monte Car,lo, data include also 
Once again we 



observe a reasonable accordance between theory and ex- 
periment. For g > 8 the comparison is hampered by un- 
certainties of the experimental parameters g and E c . The 
parallel conductance of earlier experiments was system- 
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FIG. 5: Coulomb oscillations for sample VII (g = 4.75, 
thin lines) compared with Quantum Monte Carlo calculations 
(thick lines). The conductance G normalized to the high tem- 
perature conductance G c i is shown as a function of the di- 
mensionless gate voltage n g = (C g U g )/e for the temperatures 
as given in the right margin. The calculations were done at 
temperatures (from bottom to top) 0.048, 0.1, 0.2, 0.22, 0.25, 
0.33, 0.4, 0.5, 0.67, 1.0 and 2.0 E c /k B . 
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FIG. 6: Maximum and minimum linear response conductance 
Gmin and G max for sample VII (g = 4.75) normalized to the 
high temperature conductance G c i- The experimental data 
(•) are compared to the results of our PIMC simulation (□). 
Data points for G max at the lowest temperatures have been 
omitted (see text). The inset shows the same data on a log- 
arithmic scale for better comparison of G m i n . Errorbars are 
only shown if they exceed the symbol size. 
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FIG. 7: Maximum and minimum conductance for l/{/3E c ) = 
0.05 as a function of the tunneling strength g for all examined 
samples (•) in comparison with PIMC calculations ^ (□ ) . In- 
cluded are also the experimental data of Joyez et ajJ (■) and 
the resultsi-ef perturbation theory in second orderHu (o) and 
third ordeicij (o). Errorbars are only shown if they exceed the 
symbolsize. 



atically underestimated by the assumption of symmetry 
of the SET while for the determination of the charging 
energy E c sufficient high temperature data were missing. 

Also shown are-tie predictions of ihe perturbation ex- 
pansion in secondolj and third ordera in g. The range of 
validity of the second-order perturbative approach is lim- 
ited to g < 2.5 where the maximum conductance G max 
drops with increasing g. The plateau and the following 
increase of the maximum conductance can not be de- 
scribed by perturbation theory. Also for G min deviations 
occur at g > 2.5 while the Monte Carlo approach gives 
excellent results up to g = 10. 



V. CONCLUSION 

We have presented experimental results for the con- 
ductance of single-electron transistors as a function of 
temperature and dimensionless gate voltage. The em- 
ployed four-junction layout for the SET allows for an 
unambiguous determination of the physical parameters 
g and E c . Thus we were able to clarify and eliminate 
problems encountered in earlier experiments. In partic- 
ular, comparison with theory can be made without any 
adjustable parameter. We have compared the experimen- 
tal findings with perturbation theory in second order in 
g and with results of PIMC simulations. 

Comparison with perturbation theory was made for 
the maximum and minimum conductance of five SETs 
with different tunneling strength. At g < 1.85 we found 
good agreement for the maximum linear conductance 
with second-order perturbation theory for the whole 
range of temperatures. Surprisingly even for such low- 
conductance SETs deviations from perturbation theory 
for G m i n are pronounced at low temperatures. In con- 
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trast, a detailed comparison with PIMC data for g = 4.75 
revealed good agreement between experiment and the- 
ory outside the perturbative regime. Further comparison 
showed that at g — 4.75 not only G max and G m ; n but also 
the form of the Coulomb peaks could be described very 
well by our simulations for temperatures 1/((3E C ) > 0.05. 

Finally we have presented a comparison of experimen- 
tal and theoretical results for the maximum and mini- 
mum conductance as a function of the tunneling strength 
g for a fixed temperature. The experimental data for this 
comparison are combined from experiments on different 
SETs including also the results of earlier experiments by 
Joyez et al.. For g < 8 we find that the PIMC data is in 



good accordance with experiment whereas perturbation 
theory in second order shows significant deviations for 
g > 2.5. Also the available data of third-order perturba- 
tion theory still show rather large deviations for G max . 

The good agreement found in our study is an affir- 
mation that the path integral formulation in combina- 
tion with the Monte Carlo method allows for an accu- 
rate description of the minimum conductance G m i n over 



the whole range of experimentally accessible parameters. 
Limitations of the Monte Carlo method for the calcu- 
lation of the maximum conductance G max exist at low 
temperatures and small conductances due to slow con- 
vergence. Nonetheless G max as well as the entire shape 
of the Coulomb peaks can be described well in a range of 
parameters which lies outside the perturbative regime. 
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